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Abstract 

The growth and dynamics of sohd surfaces displays a multitude of power law rela- 
tionships, which are often associated with geometric self-similarity. In many cases 
the mechanisms behind these power laws are comparatively trivial, and require lit- 
tle more than dimensional analysis for their derivation. The information of interest 
to surface physicists then resides in the prefactors. This point will be illustrated 
by recent experimental and theoretical work on the growth-induced roughening of 
thin films and step fluctuations on vicinal surfaces. The conventional distinction 
between trivial and nontrivial power laws will be critically examined in general, and 
specifically in the context of persistence of step fluctuations. 
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Empirical science is apt to cloud the sight, and, by the very knowledge of functions 
and processes, to bereave the student of the manly contemplation of the whole. The 
savant becomes unpoetic. Ralph Waldo Emerson 

Per Bak's theory of self-organized criticality is based on the observation that 
power laws are ubiquitous in nature, and the understanding that this observa- 
tion requires a (general) scientific explanation. The latter point is often taken 
for granted in statistical physics, but it is less self-evident in other disciplines. 
In fact, it could be argued that the predelection of statistical physicists for 
power law relationships is a professional deformation, which originates in the 
success story of the theory of equilibrium critical phenomena. In that context, 
power laws are indeed anomalies which require the ingenious machinery of the 
renormalization group for their explanation. 

However, not every power law carries a deep message. Statistical physicists 
account for this by distinguishing between trivial and nontrivial power laws. 
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Though deeply rooted in our jargon, this distinction hard to make precise. 
Most people would agree that the relation 
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for the mean square displacement of a random walker (really a manifestation 
of the central limit theorem) is a trivial power law, whereas, say, the Onsager 
exponents for the two-dimensional Ising model are nontrivial. But consider 
Kolmogorov's 1941 theory of fully developed turbulence [1]. The derivation of 
the A;^^/^-energy spectrum is trivial, in the sense that it uses only dimensional 
analysis, but it requires the higly nontrivial physical insight that energy dis- 
sipation is constant across scales. Clearly the demarcation line between the 
trivial and the nontrivial also evolves with the progress of science. 

Since the discovery of the fractal nature of diffusion- limited aggregation (DLA) 
more than two decades ago [2] , scaling concepts have become one of the main 
tools in the study of growth processes in condensed matter [3,4,5,6]. While 
much of the theoretical activity has been driven by the (still unfinished) quest 
of understanding the two most prominent nontrivial power laws in the field 

- the fractal dimension of DLA, and the strong coupling exponents of the 
Kardar-Parisi- Zhang (KPZ) equation [7] in dimensions larger than one [4,8,9] 

- experimentalists have begun to routinely employ the concepts of scaling and 
self-similarity to analyse topographic data in thin film and crystal growth. In 
the following I will describe some recent applications of scaling ideas in sur- 
face physics. 1 will argue that, more often than not, power laws that, by the 
standards of statistical physicists, are quite trivial, have been of most use in 
interpreting experimental data and gaining insight into kinetic processes at 
real surfaces. The discussion will include the currently popular concept of per- 
sistence of a stochastic process [10], which provides an interesting perspective 
on the distinction between trivial and nontrivial power laws. 

Our first example is the growth-induced roughening of thin films of an organic 
semiconductor, diindenoperylene (DIP) [11]. The relevant scaling law is the 
relationship between the root mean square surface roughness a (the standard 
deviation of the film height distribution) and the film thickness D, 



which defines the roughening exponent /3. The experimental data for a{D) 
shown in Fig.l are remarkable in several respects. First, they cover more than 
two decades in film thickness, and they contain consistent results obtained 
using two completely different methods, x-ray diffraction and atomic force 
microscopy. The second remarkable feature is the experimentally determined 
value P — 0.75 ± 0.05 of the roughening exponent. This value exceeds the 
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Fig. 1. Experimental data for the surface roughenss a as a function of DIP film thick- 
ness. Full data points were obtained from atomic force microscopy, and open data 
points from x-ray diffraction. Inset shows a typical fit of the measured scattering 
intensity to the form expected for a rough surface (from [11]) 

random deposition (RD) limit /3rd = 1/2, which apphes in the absence of any 
mass transport between different layers of the growing film [3,6]. In this limit 
the height follows a Poisson distribution, and (2) becomes equivalent to the 
random walk relation (1). Random deposition can be approximately realized 
in the growth of crystalline metal films, where transport between atomic layers 
is inhibited by step edge barriers (see below). Since mass transport processes 
along the growing surface are driven by differences in bonding energy, they 
will generally tend to smoothen the film. The random deposition roughness 
(Trd would therefore be expected to define an upper bound on the roughness 
that random fluctuations can induce during the growth of a thin film; in 
particular, [3 < /?rd = 1/2. In this sense systems with (3 > 1/2 are anomalous, 
and constitute examples for the (largely unexplained) phenomenon of rapid 
roughening [4,11]. 

But do the DIP data in Fig.l really violate the RD bound? To answer this 
question, it is essential to include also the prefactor of the power law rela- 
tionship (2) in the analysis. Denoting by d the thickness of a single atomic 



layer, the random deposition roughness is ctrd = d^D/d, which is depicted 
as a dotted line in Fig.l. The comparison shows that the experimental rough- 
ness data remain helow the random deposition limit even for the largest film 
thicknesses. It is conceivable (though perhaps not very likely) that the rough- 
ening exponent crosses over to a value below 1/2 before the critical thickness 
Dq in Fig.l is reached, and that the early time value (3 ~ 0.75 is merely a 
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transient ^ ; unfortunately films with thicknesses D > could not be grown 
for technical reasons. If, on the other hand, the measured value j3 ^ 0.75 is 
truly asymptotic, then the observed roughness cannot be due to the random 
fluctuations in the deposition beam. In [11] it was conjectured that quenched 
in-plane disorder arising from the boundaries of the tilt domains in the organic 
thin film could be responsible for the rapid roughening behavior. Such disorder 
is known to induce sublinear roughening behavior of the form a ^ D/ ln{D)^, 
which can mimic a power law with exponent 1/2 < /? < 1 over extended time 
scales [13]. Unfortunately, the available experimental information about the 
growth process on the molecular level is insufficient at present to validate or 
refute this hypothesis. 

If our understanding of DIP film growth, as described above, seems unsatis- 
factory, it is nevertheless fairly representative of the field as a whole. Although 
the power law relationship (2) has been, and is currently being reported for 
a host of growth systems [3,5], including metallic as well as semiconductor 
materials and crystalline as well as amorphous films, and although the mea- 
sured values of (3 tend to cluster around numbers that can be derived from 
theoretical models [14] , a clear-cut example, where a well-understood growth 
mechanism gives rise to a definite prediction for /3 which is quantitatively con- 
firmed by experiment, is so far missing. In particular, 2-|-l-dimensional KPZ 
scaling, a mathematical object under passionate theoretical pursuit for close 
to twenty years, is still to emerge in any real growth system ^ . 

Perhaps the only exception to this statement is mound formation in homoepi- 
taxial crystal growth [6], where the exponent f3 describes the roughening of 
a morphology that, rather than being scale- invariant, displays a distinct lat- 
eral length scale. The growth of "wedding cakes" under conditions of strongly 
inhibited interlayer transport is particularly simple [16,17]. In this limit the 
roughening exponent takes on the random deposition value j3 = f3^x) = 1/2, 
while deviations of the height distribution from the ideal Poisson form show 
up in the prefactor. The Poisson distribution is cut off at large heights because 
a new layer can be nucleated on top of a mound only when the top terrace 
has reached a critical size [18]. This leads to the expression 



for the surface roughness, where 6c is the coverage corresponding to the crit- 
ical top terrace size [19]. Equation (3) is our first example of a trivial power 
law where the nontrivial information (the specifics of the interlayer transport 
processes that determine 6^ [18,19]) resides in the prefactor. 

^ Similar transients have been observed in simulations of metal epitaxy [12]. 

^ In 1-|-1 dimensions KPZ-scaling with the (trivial or nontrivial?) exponent /? = 1/3 

has been demonstrated for one-dimensional slow combustion fronts in paper [15]. 
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In the remainder of the paper we wiU be concerned with the roughening of 
one-dimensional objects (hnes). Rough hnes appear naturally as atomic steps 
on vicinal crystal surfaces [20,21]. A vicinal surface is obtained by cutting the 
crystal in a direction close to (in the vicinity of) a high symmetry plane, and it 
consists of high symmetry terraces separated by steps of monoatomic height. 
In thermal equilibrium, the steps are roughened by thermal fluctuations. We 
describe the step at time t by a function y{x,t), where the x-axis is taken 
along the mean step direction and the y-direction is perpendicular to the step 
(in the direction of vicinality). Up to the length scale where coUisions with 
neighboring steps become important the static step conformation is the graph 
of a one- dimensional random walk, {{y{x,t) — y{x',t)Y) ~ |a; — x'\. Time- 
dependent fluctuations are well described by a linear Langevin equation of 
the form 

d f 

-y{x,t) = -Ki-—j y{x,t) + 7^{x,t) (4) 

where X is a positive constant, r){x,t) is white noise, and the dynamic exponent 
z depends on the dominant kinetic pathway through which the step fluctua- 
tions relax to equilibrium. The most important cases arc -2 = 2, corresponding 
to fast mass exchange between the step and the terrace {nonconserved kinet- 
ics) and z — A corresponding to mass transport only along the step {conserved 
kinetics) . It is straightforward to show from (4) that the temporal step corre- 
lations scale as 

C{t) = {iyix,s)-yix,s + t)f)^t'/^. (5) 

By analogy with (2) we can say that the roughening exponent of the step is 
/3= l/2z. 

The scaling law (5) has been observed in many experiments, on both metal 
and semiconductor surfaces [20,21]. The focus in these experiments has usually 
not been on the power laws as such, which (in the sense of our introductory 
discussion) are rather trivial manifestations of the simple linear dynamics (4) . 
The relation (5) is useful mainly as a classification tool, which serves to identify 
the dominant step relaxation processes. The nontrivial (and materials-specific) 
information lies instead in the temperature dependence of the prcfactor, which 
gives insight into the energy barriers governing the atomic processes at the step 
edge. Examples of such processes are shown in the left panel of Fig. 2. The right 
panel shows simulation data for the correlation function C{t), which illustrate 
the decrease of the prcfactor of the t^/^-law as the kink rounding barrier i?kr 
is increased. As shown by the detailed analysis in [22], the activation energy 
of the prcfactor depends linearly on £^kr, provided this quantity is larger than 
the kink energy (the energy cost for the formation of a kink), and therefore 
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Fig. 2. Left panel: Schematic of atomic processes at a step edge. The diffusion of 
single atoms along straight segments of the step is governed by an energy barrier Est ■ 
The detachment from a kink requires the additional energy .Edet) and the rounding 
of the kink the energy E]^,.. The distance between kinks is ^k- Right panel: Temporal 
step correlation function C{t) obtained from Monte Carlo simulations for different 
values of £^kr- Dashed lines are fits of the form C{t) = At^l^ (from [22]) 

a temperature-dependent measurement of C (t) can be used to experimentally 
determine E'kr- 

It would be premature to conclude, however, that the kinetics of step fluctu- 
ations can be reduced completely to "trivial" power laws like (5). Although 
the Langevin equation (4) is linear, it is still considerably more complex than 
a simple random walk, because the step is a spatially extended object. As 
a consequence, the step position y{x,t) at fixed a: is a non-Markovian pro- 
cess in t. The non-Markovian character of a stochastic process manifests itself 
clearly when considering the persistence probability, i.e. the probability that 
the process does not cross a particular value in a specified time interval. The 
computation of this quantity requires the knowledge of temporal correlations 
of arbitrary order. For a Gaussian process, such as the solution of the Langevin 
equation (4), all higher order correlation functions are, in principle, encoded 
in the two-point function (5), but in practice the calculation of the persistence 
probability for a general Gaussian process is a hard, unsolved problem [10] . 

To be concrete, we define the step persistence probability 

P{to,t) = Prob[|/(s) ^ y{to)\to <s<to + t], (6) 

where it is understood that the step is completly straight at time t = 
{y{x,0) = 0). Note that here the specified level that the step should not 
cross is set by its position at the beginning of the time interval. The analytic 
and numerical investigations in [23] have shown that the quantity (6) displays 
two distinct scaling regimes. In the transient regime t^ <^ t the persistence 
probability decays as P ~ t~^", while in the steady state regime to ^ ^ it 
decays as P ~ t"^® The two persistence exponents 6q and 6s are different. 
Moreover, while the steady state exponent is simply related to the roughening 
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Fig. 3. Left panel: Transient persistence exponent a-s a function of the roughening 
exponent (5 for the Hnear Langevin equation (4). The figure shows numerical esti- 
mates obtained from discrete solid-on-solid models (squares), direct simulation of 
the Langevin equation (triangles), simulation of an equivalent stationary Gaussian 
process (diamonds), a perturbation expansion around the Markovian case (3 = 1/2 
(bold line) and rigorous bounds derived by comparison with related Markovian pro- 
cesses (gray area) (from [23]). Right panel: Experimentally determined persistence 
probability for steps on a Si(lll)-surface which has undergone a reconstruction by 
deposition of a fraction of a monolayer of Al. Symbols show measurements at 970 
K (squares), 870 K (circles) and 770 K (triangles). The data sets have been shifted 
vertically for greater clarity; in fact the data show no systematic temperature de- 
pendence (from [28]) 

exponent (3 by 



/3, 



(7) 



the transient exponent is truly "nontrivial" , in the sense that it can be com- 
puted only approximately using a range of fairly sophisticated methods (see 
Fig.3). 

The relationship (7) provides a good illustration of the fallibility of the dis- 
tinction between trivial and nontrivial scaling exponents. On the one hand it 
would seem that ds-, being simply related to the arguably trivial exponent j3, 
should also be classified as trivial. On the other hand the relation (7), while 
well supported by simulations and plausible scaling arguments, is not rigor- 
ously established; in fact, in the context of the general problem of determining 
the persistence probability of a Gaussian process with known correlator, there 
is no reason why 9s should be simpler to compute, and hence less nontrivial, 
than 6*0. The scaling relation (7) holds because the stochastic process y{x,t) 
at fixed x, in the steady state regime, combines two invariance properties: It 
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is scale invariant as well as translationally invariant in time ^ [26] . Because it 
is based on these general principles, (7) appears to hold also for non-Gaussian 
interface fluctuations, such as those of the KPZ universality class, where not 
even the temporal two-point correlations (let alone correlations of higher or- 
der) are explicitly known [27]. 

Several groups have recently met the challenge of experimentally measuring 
the persistence probability of one- dimensional interfaces, for steps on vicinal 
surfaces [28,29,30] as well as for slow combustion fronts in paper [31]. In all 
cases the steady state persistence exponent 9s was measured, and the scaling 
relation (7) was confirmed, with /3 = 1/4 [28], f3 = l/S [29] and P = 1/3 [31], 
respectively. The transient persistence exponent (which is of greater theoretical 
interest) has so far not been experimentally accessible, because of the difficulty 
of preparing the special (fiat) initial condition that it refers to. In addition to 
the persistence probabihty defined by (6), in [28] also the probabihty for the 
step not to cross the mean step position (determined as an average over the 
measured time series) was considered. This quantity turns out not to follow 
a power law. This is because in the thermodynamic limit the step position 
starts, with probability one, infinitely far away from the mean (the variance 
of y diverges with diverging step length); hence the persistence probability 
remains at unity and decays exponentially only on a much larger time scale 
set by the finite length of the step [32,33]. 

Experimental step persistence data from [28] are displayed in the right panel 
of Fig. 3. In contrast to the step correlation function in Fig. 2, the persistence 
data show no systematic temperature dependence; the prefactor of the ob- 
served power law depends only on the sampling time [34]. This is because 
the property of a height fluctuation to return (or not) to its initial value in 
a prescribed time interval does not depend on the overall amplitude of the 
fluctuation; the persistence probability is a functional only of the shape of 
the (suitably normalized) correlation function [16]. As a consequence, per- 
sistence measurements cannot be used to extract energy barriers for specific 
microscopic processes. The benefit of these measurements is of a more funda- 
mental nature - they prove that the theoretical description of step fluctuations 
through the Langevin equation (4) , and the underlying picture of "universality 
classes" encoded by the dynamic exponent z, extends beyond the two-point 
function to correlations of arbitrary order. As the persistence probability is 
known to be extremely sensitive to hidden temporal correlations affecting the 
interface fluctuations [27], it may be particularly useful in cases where the as- 
signment of the universality class (the value of z) through more conventional 



^ Together these two properties define a Gaussian process known as fractional Brow- 
nian motion (fBm) [24], a generalization of tlie Wiener process, the continuum limit 
of the simple random walk. For another recent application of fBm to a persistence 
problem sec [25]. 
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measurements is ambiguous. 

Wiiere does this leave us with regard to the general reflections on power laws, 
the trivial and the nontrivial, the deep and the shallow, which introduced this 
paper? Prom discussions in the early 1990's, when SOC was challenged by 
the competing concept of Generic Scale Invariance [36], I recall that Per Bak 
had a distinct, and rather unfavorable opinion of simple power laws generated 
by simple equations such as (4): He referred to them as systems operating 
by the garbage in, garbage out principle, because they merely perform a (lin- 
ear, or, more generally, nonlinear) transformation of the driving white noise 
into correlated fluctuations. Per viewed self-organization as an essential part 
of SOC. He emphasized the capability of self-organizing systems to develop 
new, emergent levels of structure, and to undergo a history which is open to 
contingent influences. Is a rough surface or a fractal cluster grown by DLA 
self-organized in this sense? Probably not ^ . Still, we must accept that many 
of the power laws in the world that surrounds us have rather humble origins. 
As theoretical physicists, our task it to explain these origins; but as "experi- 
mental philosophers" (Per Bak) we would also like to know what they mean. 
The latter question is of course not one to be answered, but one that is to be 
constantly clarified (and reobscured) in the ongoing discourse of the commu- 
nity. In these discussions Per will be sorely missed, and gratefully remembered 
for a long time to come. 



Acknowledgements 

I wish to thank Satya Majumdar for helpful correspondence, and Arndt Diirr 
and Ellen Williams for the permission to reproduce experimental data. The 
work reported here was conducted at the University of Duisburg-Essen with 
partial support of DFG within SFB 237. 



References 

[1] U. Prisch: Turbulence (Cambridge University Press, 1995) 

[2] T.A. Witten, L.M. Sander: Phys. Rev. Lett. 47, 1400 (1981) 

^ Nevertheless it is true that scale-invariant structures need sufficient time to de- 
velop; more precisely, the time T required to grow a structure of size R typically 
scales as T with z > 1. This is probably the reason why most spatial power 

laws in Nature extend over less than two orders of magnitude [35]: Under typi- 
cal physical conditions, there is insufficient time for scale-invariant correlations to 
develop further. 



9 



[3] A.-L. Barabasi, H.E. Stanley: Fractal Concepts in Surface Growth, (Cambridge 
University Press, 1995) 

[4] J. Krug: Adv. Phys. 46, 139 (1997) 

[5] P. Meakin: Fractals, Scaling and Growth far from Equilibrium (Cambridge 
University Press, 1998) 

[6] T. Michely and J. Krug, Islands, Mounds and Atoms. Patterns and Processes 
in Crystal Growth Far from Equilibrium (Springer, Berlin 2003) 

[7] M. Kardar, G. Parisi, Y.-Z. Zhang: Phys. Rev. Lett. 56, 889 (1986) 

[8] J. Krug and H. Spohn, in: Solids far from Equilibrium, ed. by C. Godreche 
(Cambridge University Press, Cambridge 1991) pp. 479-582 

[9] T. Halpin-Hcaly, Y.C. Zhang: Phys. Rep. 254, 215 (1995) 

[10] S.N. Majumdar: Curr. Sci. 77, 370 (1999) 

[11] A.C. Diirr, F. Schreiber, K.A. Ritley, V. Kruppa, J. Krug, H. Dosch, B. Struth: 
Phys. Rev. Lett. 90, 016104 (2003) 

[12] K.J. Caspersen, A.R. Layson, C.R. Stoldt, V. Fournee, P.A. Thiel, J.W. Evans: 
Phys. Rev. B 65, 193407 (2002) 

[13] J. Krug: Phys. Rev. Lett. 75, 1795 (1995) 

[14] J. Krim, G. Palasantzas: Int. J. mod. Phys. B 9, 599 (1995) 

[15] J. Maunuksela, M. Myllys, O.-P. Kahkonen, J. Timonen, N. Provatas, M. J. 
Alava, T. Ala-Nissila: Phys. Rev. Lett. 79, 1515 (1997) 

[16] J. Krug: J. Stat. Phys. 87, 505 (1997) 

[17] M. Kalff, P. Smilauer, G. Comsa, T. Michely: Surf. Sci. 426, L447 (1999) 

[18] J. Krug, P. Politi, T. Michely: Phys. Rev. B 61, 14037 (2000) 

[19] J. Krug, P. Kuhn, in: Atomistic Aspects of Epitaxial Growth, ed. by M. Kotrla, 
N.L Papanicolaou, D.D. Vvedensky, L.T. Wille (Kluwer, Dordrecht 2002) 
pp. 145-163 

[20] H.-C. Jeong, E.D. Williams: Surf. Sci. Rep. 34, 171 (1999) 

[21] M. Giesen: Prog. Surf. Sci. 68, 1 (2001) 

[22] J. Kallunki, J. Krug: Surf. Sci. Lett. 523, L53 (2003) 

[23] J. Krug, H. Kallabis, S.N. Majumdar, S. Cornell, A.J. Bray, C. Sire: Phys. Rev. 
E 56, 2792 (1997) 

[24] B.B. Mandelbrot, J.W. van Ness: SIAM Rev. 10, 422 (1968) 

[25] S.N. Majumdar: Phys. Rev. E 68, 050101 (R) (2003) 



10 



[26] J. Krug: Markov Proc. Rel. Fields 4, 509 (1998) 
[27] H. Kallabis, J. Krug: Europhys. Lett. 45, 20 (1999) 

[28] D.B. Dougherty, 1. Lyubinetsky, E.D. Williams, M. Constantin, C. Dasgupta, 
S. Das Sarma: Phys. Rev. Lett. 89, 136102 (2002) 

[29] D.B. Dougherty, O. Bondarchuk, M. Degawa, E.D. Williams: Surf. Sci. Lett. 
527, L213 (2003) 

[30] M. Constantin, S. Das Sarma, C. Dasgupta, O. Bondarchuk, D.B. Dougherty, 
E.D. Williams: Phys. Rev. Lett. 91, 086103 (2003) 

[31] J. Merikoski, J. Maunuksela, M. Myllys, J. Timonen: Phys. Rev. Lett. 90, 
024501 (2003) 

[32] C. Dasgupta, M. Constantin, S. Das Sarma, S.N. Majumdar: Preprint 
( cond-iiiat/0 307086) 

[33] S.N. Majumdar (private communication) 

[34] E.D. Williams (private communication) 

[35] D.Avnir, O. Biham, D. Lidar, O. Malcai: Science 279, 39 (1998) 

[36] G. Grinstein, D.-H. Lee, S. Sachdev: Phys. Rev. Lett. 64, 1927 (1990) 



11 



